A simple statistical explanation for the localization of energy in nonlinear lattices 

with two conserved quantities 
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The localization of energy in the discrete nonlinear Schrodinger equation is explained with sta- 
tistical methods. The partition function and the entropy of the system are computed for low- 
amplitude initial conditions. Detailed predictions for the long-time solution are derived. Localized 
high- amplitude excitations absorb a surplus of energy when they emerge as a by-product of the 
production of entropy in the small fluctuations. The thermodynamic interpretation of this process 
applies to many dynamical systems with two conserved quantities. 
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Localization of energy within a small number of iso- 
lated high-amplitude structures is a widespread phe- 
nomenon in nonlinear opticsQ, 0, Q , plasma physics , 
Bose-Einstein condensates |5j and nonlinear lattice dy- 
namics 0, S E3 ■ Peaks of the energy density result 
either from a collapsing wave train that leads to a finite- 
time blow-up of the amplitude 0, |n|, or, in spatially 
discrete systems, from a sequence of merging breathers 

It is the purpose of this article is to show that the for- 
mation of peaks is driven by the production of entropy 
and that this behavior is generic for the thermalization 
of many conservative systems where a second quantity 
is conserved in addition to the Hamiltonian. The dis- 
crete nonlinear Schrodinger equation is a simple generic 
equation that describes discrete breathers in nonlinear 
optical waveguide arrays 0, Q and dilute Bose-Einstein- 
condensates that are trapped in periodic potentials 
The spatial discreteness avoids the leakage of energy 
to infinitesimal scales that can occur during the wave- 
collapse of continuous systems @ so that this system is 
a simple but representative model for localization pro- 
cesses. The non-compactness of the phase space leads to 
a merely technical difficulty in the statistical treatment 
of high-amplitude structures. 

Fig. la shows the focusing process for the focusing dis- 
crete nonlinear Schrodinger equation (DNLS) 



i4>n = <Pn+l + <t>n-l + \4>n\ 2 4>r, 



(1) 



Any coefficients of this equation can be removed by scal- 
ing, and consequently all quantities are dimensionless. 
Firstly, breathers with moderate amplitudes appear pe- 
riodically in space and time following a phase instability 
of a regular low-amplitude initial solution. Subsequently, 
they merge into more persistent peaks with higher am- 
plitudes (lattice site 188 in Fig. la). The system finally 
settles into a state where immobile high-amplitude peaks 
(the ring with \(f>\ w 2.3 in Fig. lb) emerge from a low- 
amplitude disordered background (core with \(j)\ < 0.5). 
The peaks divide the system into patches of the order 



of 100 lattice sites where the amplitude is low and the 
dynamics is irregular. The peaks oscillate corresponding 
to their amplitude-dependent frequency and their ampli- 
tudes fluctuate slightly, but their position in the lattice 
almost never changes. 

This behavior depends crucially on the system's energy 
E =< H >, where H = EMti+i + #&fi) + l<t>Uf 
is the Hamiltonian of the DNLS i<p n — Jp-. Persis- 
tent localization of energy occurs only if the system's 
energy is positive. In this case the peaks finally absorb 
almost the total energy (Fig. 2a). The height of the peaks 
\(j)\ w 2. ..2. 5 is almost independent of the total energy 
(Fig. 2b), but the number of peaks increases with the 
energy. For negative energies, there is no localization 
of energy (Fig.lc) and the system settles into a state of 
low-amplitude fluctuations (Fig. Id, Fig. 2a). The system's 
second conserved quantity, the modulus-square norm (or 
'particle number') A =< A > with A — 'fii'fii ^ s a ^ so 
crucial for this phenomenon. There is no persistent lo- 
calization of energy in this system if the rotational sym- 
metry linked to this second conserved quantity is broken. 
The discrete nonlinear Schrodinger equation admits en- 
ergies — 2 A < E < 2 A for low-amplitude initial condi- 
tions A/N < 1. As suggested by Fig. lb, the phase space 
may be separated into a low-amplitude domain A< for 
|(/><| < r and a high-amplitude domain A> for |0>| > r. 
To avoid formal difficulties with diverging terms, a pre- 
liminary upper bound R > r of the phase space is in- 
troduced as \4>\ < R. The main contributions to the 
energy arise from interactions in the low-amplitude do- 
main (where N—K oscillators are gathered and nonlinear 
contributions are negligible) and from the quartic contri- 
bution of K <C N — K oscillators in the high-amplitude 
domain. The Hamiltonian is therefore approximated as 

Ti. ~ 7i< + 7i> 

= X>( r - l<^l) u O- l<fc+il)(&0i+i + 

+£u(|&|-r)## a 

(2) 

with the unit step function u{x < 0) = 0, u{x > 0) = 1. 
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The second integral of motion may be separated in a 
similar way as A — A< + Ay . The system's entropy as 
a function of the conserved quantities E and A can be 
computed from the grand partition function y(/3, 7) = 
/ PYliLi d <t>i d <t>i usin g the density p = e ~m-iA). The 
partition function is a sum over all possible numbers of 
peaks K, and only its leading term 
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y< ((3, % N-K) y> (p,"/,K) (3) 
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will be considered. K and R will be determined later 
in order to maximize the entropy. The factor (^) gives 
the number of combinations of K high-amplitude sites 
on N lattice sites. The contributions y < and ?/> will 
be computed separately. The patches of small-amplitude 
fluctuations between the peaks give the contribution 
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FIG. 2: Energy (dots in (a)) and height (dots in (b)) of peaks 
with \<j>\ > 1 as a function of the total energy after 2 x 10 time 
steps of numerical integration. The initial conditions of the 
4096 lattice sites are waves with an amplitude \<f>\ = 0.3 and 
wavenumbers from k = to k = it. Lines: Thermodynamic 
equilibrium value (a) and critical height for lattice pinning 
(b). 



(4) 



to the partition function with p < = e~ l3 ( H< ~~ lA< \ These 
regions are in thermal equilibrium with each other, and 
therefore they are regarded as one single held of N — K 
oscillators. Similarly, the K oscillators in the high- 
amplitude regime A> have degrees of freedom corre- 
sponding to variations of phase and amplitude of the 



peaks. This leads 
function 



y>(/?,7) 



to the contribution to the partition 
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(d) 



term of the peaks which is small compared to the quar- 
tic energy. However, peaks and fluctuations are coupled 
thermally, so that they can exchange energy and parti- 
cles. This exchange vanishes on average when peaks and 
fluctuations have the same temperature /3 _1 and chem- 
* ical potential 7. The energy E w E < + E> and the 
particle number A = A < + Ay have contributions from 

* / either domain. A, A <7 Ay, E> are positive, E and -E< 

t ' may be positive or negative. 

The partition function y < can be reduced to Gaus- 
Ri^rn integrals. The density may be written as 
^ = rj e -( Al ^- A2 ^+i) 2 e-( Al ^'- A2 ^'+i) 2 with A1/2 = 
(^-/3(7 + 2) ± y/-/3(j- 2))/2 (it turns out later that 
—Pi > and I7I > 2). Introducing the variables 
x n = Xx4>' n — ^20^+1, the partition function of the flue- 



tuations is reduced to Gaussian integrals 

V<(P,l) = {( e~< Aet{d^Jdx n )\{dx n ) 2 



(6) 



Using det(dx n /d(j)' n ) = A 



N-K 



± A 



N-K 



the Jacobian is 



-2.C- 



_Mk(d4>'Jdx n ) w Ai {N K) for N - if > 1 since Ai > A 2 . 
iTne square in equation @ is obtained from the identical 
^integration over tfi" . This gives an analytic expression for 
the partition functio n y < (P,j, N -K) = A N ~ K with A = 
tt\^ 2 = 7r(-/37 + V/? 2 (7 2 _ 4)) _1 . The corresponding 
energy is 
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FIG. 1: Numerical integration of the DNLS with 4096 oscil- 
lators. The initial conditions are waves with the wavenumber 
k = (a),(b), and with k = tt/2 (c),(d) for \4>\ = 0.3. (a),(c): 
Spatiotemporal pattern of high-amplitude states (dark gray) 
in a small sector of the chain for the first 2000 time steps. 
(b),(d): Distribution of cj> after 2 x 10 s time steps. 
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(7) 
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and their particle number is 
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(8) 



Consequently, the inverse temperature is f3 = —2E < (N — 
i4f)/(4A^ — E 2 ,) and the chemical potential is 7 = 
{AA\ + E\)/{2E < A < ). Using these expressions, the 
canonic transformation S < — lny < + j3{E < — 7^<) leads 
to the entropy of the small fluctuations 



(N-K)lntt 



(9) 



with n = {4A 2 , - Ei-)/(4A < (N - K)). The entropy 
is the most useful thermodynamic potential in this con- 
text because its variables (particle number and energy) 
are known from the initial conditions. The fluctuations 
entropy per lattice site Infl is plotted in Fig. 3a. The 
lines E < = 2A < and E < = —2A < correspond to waves 
with k — and k — n respectively. The entropy is in- 
finitely low for these ordered states. The ridge E < = 
corresponds to a fluctuating state with an infinite tem- 
perature where all wavenumbers have the same power. 
A wave with k = n/2 is a non-thermalized solution with 
E = 0. This line (3 = was first computed in [jj where 
the additional nonlinear correction were included, and it 
was identified as the transition line to the localization 
phase. No analytic results for the statistics beyond this 
line have been available yet. 

Equation gives a valid expression for the total en- 
tropy if the system's total energy is negative (corre- 
sponding to the right slope with E < < of Fig. 3a). 
The temperature is positive in this case, and conse- 
quently the density p> ~ e ~' 3 ^-'l^"l I 2 decays rapidly 
for huge amplitudes so that the high-amplitude contribu- 
tion to the partition function is negligible. In this phase 
only small fluctuations contribute to the total entropy 
S = Nlnfl. This describes the low-amplitude fluctuating 
state (E = E K ,A = A<) with no peaks of Fig. Id. 
The slope E < > is linked to negative temperatures 
0- The density p> increases with the amplitude in this 
regime which was suggested to be the reason for the for- 
mation of high-amplitude structures Q. The increase 
of the density causes obvious technical difficulties: The 
phase space is (unlike in the spin system of H,E3) non- 
compact, so that the integral y> will diverge for any fixed 
negative values of (3. However, the artificial upper bound- 
ary R of the phase space prevents this divergence. The 
partition function can first be computed as a function of 
this parameter R, which subsequently will be allowed to 
go to infinity as this maximizes the entropy. From this 
the contributions of the peaks to the entropy, the energy 
and the particle number can be computed. 
The grand partition function of the peaks increases as 
V>{P,i) ~ (e-^ 4 / 2 -^ 2 ' /(-[3R 2 + 0j)) K with the cut- 
off R of the phase space if R is greater than y/y. Since the 




0.1 



0.15 

(K/N)ln r 



E</N 
(b) 




0.05 



A</N 



E<//V 



FIG. 3: Inil (a) and (K/N)lnT (b) as a function of E < and 
A < ; line R 2 = 7 for the simulation of Fig.la,b 



density p gathers at the border R, the particle number 
related to the domain A> is Ay = KR 2 and the energy 
is E y = KR A /2. Inserting 0, 7 and R = X /2E > /A > into 
the entropy of the peaks 5> ~ —Kln(—(3R 2 ) yields 



Si 



Kln( f ( 4A <-^) l 
Kln(K/N) + Klnn 



KlnT 



(10) 



with T = 2A < N/A > E < . The combinations of K = 
A 2 > /2E > peaks on N lattice sites yield an entropy con- 
tribution S p = Iti(k) ~ Kln{N / K). Adding up S < , 5>, 
and S p , the total entropy is 



S = Nlnfl 



A 2 

A >-lnT 



2E-, 



(11) 



This expression describes the thermodynamics in the 
domain where localization occurs. The small fluctu- 
ations provide the leading term Nln£l of the entropy 
while the contribution KlnT (Fig. 3b) from the peaks is 
negligible under the constraint R 2 > 7 or 2E > /A > > 
(4A 2 . + E 2 < )/(2E < A < ). y > and V can be computed in 
the same way for various types of nonlinearities. The 
entropy increases with A < and decreases with |i?<| and 
has its maximum at E < = 0, A < = A, E > = E, Ay = 0. 
The state of maximum entropy is related to low ampli- 
tude fluctuations with /3 = and ^7 = — N/A. 
The crucial point is, that the high peaks contribute little 
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FIG. 4: Isentropes Si > S 2 > S 3 > ... for the left slope 
of Fig. 3a; trace (E < ,A < ) of the integration of Fig.la,b over 
£'<, A<; trace for a version of the DNLS where the second 
conservation law is violated. 



to the total entropy, while they can absorb high amounts 
of energy using only few particles. On the other side, the 
fluctuations can reach a state with a maximum entropy, 
if they contain the ideal amount of energy. This shows 
the thermodynamical nature of energy-localization: In 
order to maximize the system's total entropy, the ideal 
amount of energy E < must be allocated to the fluctua- 
tions. Starting from an initial condition with a positive 
energy E = E < at the left slope in Fig. 3a, the entropy 
can be increased when E < decreases while the released 
energy is stored in the localized structures. The state of 
maximum entropy corresponds to only one peak which 
absorbs the total energy (Fig. 2a) while consuming very 
few particles. This also shows the self-consistency of the 
truncation J5J that neglects the interaction of the peaks 
with their environment. The equilibrium state with local- 
ized structures corresponds to (3 = 0, so the temperature 
is not negative. In the opposite case with negative total 
energies E and a positive temperature, there is no energy 
surplus and consequently no localization (Fig.lc,d). 
Fig. 4 gives a numerical picture of the shift of energy 
and particles into the peaks for equation Q and the 
corresponding picture for a version of the DNLS where 
the rotational symmetry is broken by a small term 
O.O15i?e(0„). It shows the evolution of the fluctuations 
energy E < and particle number A < versus the corre- 
sponding equilibrium isentropes S = const for a homoge- 
neous low-amplitude (\<fi\ — 0.3) initial state E < = E = 
2A and A < = A. The boundary between A< and A> is 
r = 1. The particlc-nonconserving system produces no 
persistent peaks. Instead, it increases its entropy by in- 
creasing the total number of particles and ends up in 
a fluctuating high-amplitude state 0]. The particle- 
conserving system approaches the entropy-maximum 
by generating peaks and transferring energy from the 



fluctuations E< to E > . On the other side, growing peaks 
also absorb particles that must be shifted from A< to 
Ay. The trace (E <1 A < ) therefore strives down to the 
left in Fig. 4. This loss of particles in the fluctuations is 
unfavorable for the production of entropy. The merging 
of peaks however allows the system to store more energy 
in the peaks while feeding particles from the peaks back 
into the fluctuations A < which leads to an additional in- 
crease of the entropy. 

A state with only one huge peak is not reached experi- 
mentally (Fig. lb). Instead, a number of coherent struc- 
tures of moderate height (Fig. 2b) survive even on very 
long time scales (10 7 time steps in numerical simula- 
tions). In Fig. 4, the trace (E <1 A < ) does not reach the 
state of maximum entropy at (E < — 0, A < w 370), but 
ends up at in a state with only half that particle con- 
tent, and still a positive amount of energy. The coherent 
structures contain a significant amount of particles, and 
the entropy of the fluctuations is below its maximum. 
The reason for this is the eventual breakdown of the two 
entropy-enhancing mechanisms: Firstly, the merging of 
peaks becomes impossible because peaks above a certain 
critical height are immobilized by a lattice-pinning ef- 
fect. For that reason the trace (E < ,A < ) in Fig. 4 does 
not cross the isentropes S3, S2, S\. Secondly, the growth 
of peaks stops when this leads to no further increase of 
the fluctuations entropy. This happens when the entropy 
gain due to the energy transfer to the peaks is matched 
by the entropy loss due to the particle transfer. For this 
reason, the trace (E < ,A < ) in Fig. 4 does not approach 
the line E < — any further after ~ 10 5 time steps. 
The mobility of a peak with an amplitude \<f> n \ = \fa at 
a site n depends on the compatibility of the local con- 
servation of its energy E = a 2 /2 and its particle number 
A = a during a possible migration to the site n + 1 . Most 
of the a particles are gathered at these two lattice sites 
during the migration. The trajectory of such a migra- 
tion is therefore close to the intersection of the level set 
of E and A with the constraint of nonvanishing ampli- 
tudes at n and n + 1 only. These paths are given by 
<t> n (v) = VWT^e 1 * and <f> n+ i(v) = y/aue^e** with 
cos a = a\/ v — v 2 /2. The parameter v goes from to 1 
during the migration as the particles are shifted from n to 
n+1. ip is a phase factor. The bottleneck of this process 
is the intermediate storage of |0| 4 -energy in the coupling 
term 2i?e(^„0* +1 ): The solvability condition cos a < 1 
for all v requires that a < 4, so that this migrational 
path only exists for peaks below a maximum amplitude 
\4>n{v = 0)| < 2. For higher peaks, it is impossible to 
conserve both particle number and energy at the instant 
when \4> n \ = |<^n-i-i|- These idealized migrational paths 
deviate from exact solutions of the DNLS which necessar- 
ily have nonvanishing amplitudes for adjacent oscillators 
at n — 1, n + 2, etc. Monte-Carlo simulations of paths 
that include small amplitudes for adjacent oscillators (4- 
20 lattice sites) show that migrations for slightly higher 
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peaks {\4> n {v = 0)| < 2.28, Fig. 2b) are permitted by the 
conservation laws. Higher peaks first need to decrease by 
transferring particles to remote lattice sites before they 
can move. This defocusing process would require a very 
unlikely decrease of entropy and the local conservation 
of particles and energy is therefore statistically favor- 
able. Spontaneous collapses of peaks are only possible 
if this increases the total entropy, which is the case if the 
fluctuations have small wavelengths and a positive tem- 
perature. 

Further growth of the peaks requires a flow both of 
particles and of energy from the fluctuations to the 
peaks as the number of pinned peaks is fixed on rel- 
evant time scales. This is thermodynamically favor- 
able only if the increase of the entropy due to the en- 
ergy transfer is bigger than the decrease caused by the 
particle transfer. This process stops when the trace 
(E < ,A < ) in Fig. 4 approaches an isentrope tangentially as 
dE > /dA > \ K =conat = dE < /dA < \ s =const which is equiv- 
alent to 7 — R 2 . On this line, growth or decay pro- 
cesses of peaks absorb or release energy and particles in 
a ratio that amounts to isentropic changes of the fluc- 
tuations. The statistical results reflect microscopic dy- 
namical findings 0] of the growth and decay processes 
of localized structures perturbed by one or two incoming 
waves. Long waves lead to growth, short waves to de- 
cay of a peak while the radiated harmonics increase the 
systems entropy. Experimentally, one still finds irregular 
oscillations of the peak heights, but no average growth. 
Interactions of peaks and fluctuations that increase the 
peaks are matched on average by those interactions that 
decrease the peaks. Waves with all wavenumbers coexist 
in the state of equilibrium where growth processes are 
balanced by decay processes. Any changes of the peak 
amplitude are statistically unfavorable as they decrease 
the total entropy. 

To conclude, localization in non-integrable systems con- 
strained by two integrals of motion is a statistical process. 
The entropy is dominated by contributions from small- 
amplitude fluctuations. The entropy is maximal if an 
optimal share of each conserved quantity is allocated to 
the fluctuations. There is no localization if not enough 
energy is supplied by the initial conditions. If a surplus of 
energy is provided by the initial conditions, it is dumped 
into high-amplitude structures that absorb high amounts 
of energy (Fig. 2) while using few particles. This explains 
the phase where energy-localization occurs. The lattice 
pinning effect related to the conservation laws prevents 
the system from reaching the absolute entropy maximum. 
For the final peak size, growth- and decay-interactions of 
the peaks with the fluctuations are balanced. 
Obviously, the statistical analysis shows the macroscopic 
properties of an ensemble of microstates, and one might 
expect that the Arnold diffusion process transfers the tra- 
jectory from any initial condition to this most probable 
state. However, the thermalization may have extremely 



long transients. In the numerical simulations, the ampli- 
tudes in the initial conditions were small enough so that 
the separation of the partition function in the statistical 
description was valid, but large enough for the nonin- 
tegrability to have an impact on moderate time scales. 
In the 'integrable' limit of smooth initial conditions with 
very low amplitudes, the dynamics is most similar to the 
continuous oncdimcnsional nonlinear Schrodinger equa- 
tion, which is integrable. Consequently, this part of the 
phase space is partitioned by Kolmogorov-Arnold-Moser 
tori that are not destroyed by the nonintegrability and 
the dynamics is quasiperiodic on very long time scales. 
To escape from this state, the amplitude has to reach a 
sufficient height so that the nonlinearity absorbs a sub- 
stantial amount of energy, which becomes a rare event for 
small average amplitudes. Such quantities that are al- 
most conserved on moderate time scales can also be rele- 
vant for localization effects in systems where the Hamilto- 
nian is the only exactly conserved quantity. For instance, 
the DNLS with broken rotational symmetry shows en- 
ergy localization on shorter time scales where the particle 
number changes very little. This is the case in the early 
stadium of the path E < , A < where A is not conserved 
in Fig. 4. It is an interesting question if this also applies 
to other breather systems where the Hamiltonian is the 
only conserved quantity 0, Q] . 

The mechanism of localization is found in very diverse dy- 
namical systems as it relies only on the properties of the 
entropy functional and on the existence of two conserved 
quantities and not on the spatial discreteness. Depending 
on the type of the nonlinearity, the localized structures 
may absorb primarily the second conserved quantity Q 
and not the energy. In continuous systems described by 
partial differential equations , the transport of fluc- 

tuations is continued down to the molecular scale. Again, 
the the conservation laws require the formation of local- 
ized structures for this exploitation of degrees of freedom 
on short space scales. 
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